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(g) Error detection in seismic data 

(57) The present invention relates to the detection 
of errors in seismic data and in particular to 
apparatus for detecting errors in frames of seis- 
mic data that can be presented as images. Due 
to the techniques used for collecting seismic 
data, an image produced to represent that date 
will typically be expected to exhibit an inherent 
symmetry about an axis. The apparatus of the 
present invention comprises a storage device 
for storing a frame of seismic data which can be 
presented as an image, a selection device for 
selecting an axis within the image as the axis of 
symmetry, and a divider for dividing the date 
about the axis of symmetry to create first and 
second sets of data. Further it includes a com- 
parison means for comparing the sets of data 
according to preselected matching criteria to 
determine, based on the degree of matching, 
whether the extent of mismatch in the frame of 
data is within a predetermined tolerance. 
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The present invention relates to the detection of 
errors in seismic data and in particular to apparatus 
for detecting errors in frames of seismic data that can 
be presented as images. Due to the techniques used 
for collecting seismic data, an image produced to rep- 5 
resent that data will typically be expected to exhibit an 
inherent symmetry about an axis. 

However seismic survey data is subject to vari- 
ous errors from many causes; the most frequent 
cause of errors are background noise and misplaced 10 
equipment Analysis of the seismic data is typically a 
very long process with the data passing through many 
stages of computational processing. Hence it is desir- 
able to identify errors in frames of data as early as 
possible to avoid wasting valuable time processing 15 
bad data. A reliable technique for detecting errors 
could save much re-work. 

Typically a company involved in obtaining and 
analysing seismic data will have to transport a great 
deal of hardware to the survey location, will carry out 20 
a large number of seismic tests, and then will return 
to the laboratory to analyse the data. Seismic tests 
usually involve positioning cables in lines on the sur- 
face of the land (or towing cables behind a ship on the 
surface of water), each cable having a number of sen- 25 
sors attached to it along its length. A shock wave is 
then generated, eg. by firing a shot, at a location on 
the surface symmetrically disposed with respect to 
the cables, and the shock waves are detected along 
these cables of sensors. Given this type of arrange- 30 
ment, it is apparent that viewing the collections of sig- 
nals as an image and checking for symmetry is a 
good way of detecting errors. 

However, up until now, this has been done by ge- 
ophysicists checking by eye. Obviously this is very 35 
time consuming, labour intensive, and prone to er- 
rors, and so it would be highly desirable to reduce the 
time that a geophysicist must spend inspecting the 
data. The seismic industry has applied theories of ge- 
ophysics in an attempt to solve the problem of error 40 
detection, but has been unable to solve the problem 
purely within the domain of geophysics. 

It is hence an object of the present invention to 
provide an improved technique for detecting errors in 
frames of seismic data, which requires much less hu- 45 
man activity. 

Ideally it would be useful if the detection of errors 
could be carried out at the survey site so that if nec- 
essary, data can be re-acquired before the testing 
equipment is removed from the site. so 

Accordingly the present invention provides an ap- 
paratus for detecting errors in frames of seismic data 
that can be presented as images, such images being 
expected to exhibit an inherent symmetry about an 
axis, the apparatus comprising: a storage device for ss 
storing a frame of seismic data which can be present- 
ed as an image; a selection device for selecting an 
axis within the image as the axis of symmetry; a div- 



ider for dividing the data about the axis of symmetry 
to create first and second sets of data; and a compar- 
ison means for comparing the sets of data according 
to preselected matching criteria to determine, based 
on the degree of matching, whether the extent of mis- 
match in the frame of data is within a predetermined 
tolerance. 

Typically the comparison means will comprise a 
matching device for comparing the sets of data to de- 
termine the degree of matching between such sets, 
and decision logic, responsive to output from the 
matching device, to determine whether the extent of 
mismatch is within the predetermined tolerance. 

In preferred embodiments the comparison 
means will also comprise an edge detector for deter- 
mining the data points representing discontinuities in 
the image, and to create an edge list of such data 
points. Such edge lists will then be used as the basis 
of comparison by the matching device. The advan- 
tage of using an edge list as the input to the matching 
device is that it will typically present many less points 
for comparison than the original data would present 
(eg 2058 edgels in a 128 x 128 image (16384 bytes)). 
Hence the speed of the comparison process is im- 
proved. 

The divider can be arranged to operate on the 
whole frame of seismic data prior to passing the first 
and second sets of data separately through the edge 
detector. In this arrangement the edge lists created by 
the edge detector will be received by the matching de- 
vice for comparison. However in preferred embodi- 
ments the edge detector is arranged to operate on the 
whole frame of seismic data, and then the edge list is 
operated on by the divider so as to create two edge 
lists as the first and second sets of data for compar- 
ison by the matching device. This latter approach re- 
duces computing overhead since only one set of data 
has to be passed through the edge detector rather 
than two; the edge detection process has a signifi- 
cant setup overhead. 

The divider may just divide the data about the axis 
of symmetry, or additionally may include a mirroring 
means for re-ordering the data in the first set of data 
so that it can be represented by an image which would 
be the mirror image of the image previously repre- 
senting that set In preferred embodiments the divider 
includes a mirroring means and the comparison 
means is then arranged to compare the second set of 
data with the re-ordered first set of data. By employ- 
ing a mirroring means in the divider, a much simpler 
form of comparison means can be used than would 
otherwise be necessary. 

Typically the mirroring means re-orders the first 
set of data created by the divider by employing a buf- 
fer in which to store data from that set of data consti- 
tuting a row of an image, the contents in the buffer 
then being read out in reverse order and stored as the 
re-ordered first set of data. 
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The selection device can be implemented by 
means of an amplitude peak detector which selects 
the axis of symmetry as the position in the data hav- 
ing the highest detected intensity. However in prefer- 
red embodiments the selection device uses informa- 5 
tion contained in the seismic data, giving the position- 
ing of the shot for that frame of data, to select an axis 
as the axis of symmetry. 

The edge detector can use any one of a number 
of edge detection techniques. However in preferred 10 
embodiments the edge detector employs a Canny 
edge detection algorithm in order to identify the data 
points representing edges. The Canny algorithm is a 
known edge detection algorithm and so will not be de- 
scribed in detail herein; the article by J F Canny enti- is 
tied "A Computational Approach to Edge Detection", 
IEEE Trans. Pattern Analysis & Machine Intelligence, 
6:679-698 (1986), describes the algorithm in detail. 

In preferred embodiments the matching device 
uses a stereo matching algorithm. A particularly pre- 20 
ferable stereo matching algorithm is the PMF algo- 
rithm, which has shown in the present circumstances 
to be more efficient and robust than alternative tech- 
niques. 

Typically the seismic data represents geological 25 
structure created from seismic disturbance measure- 
ments taken at fixed points from a source of distur- 
bance. 

Viewed from a second aspect the present inven- 
tion provides a method of detecting errors in frames 30 
of seismic data that can be presented as images, 
such images being expected to exhibit an inherent 
symmetry about an axis, the method comprising the 
steps of: (a) storing a frame of seismic data that can 
be presented as an image in a storage device; (b) se- 35 
lecting an axis within the image as the axis of symme- 
try; (c) employing a divider to divide the image about 
the axis of symmetry to create first and second sets 
of data; and (d) comparing, in a comparison means, 
the sets of data according to preselected matching 40 
criteria to determine, based on the degree of match- 
ing, whether the extent of mismatch in the frame of 
data is within a predetermined tolerance. 

The present invention solves the error detection 
problem by utilising a number of techniques devel- 45 
oped for image processing to create an apparatus 
which can quickly and reliably identify errors with the 
minimum of user interaction. The apparatus 
will , fiag > frames of seismic data which show devia- 
tions from symmetry, after allowing for certain natu- so 
ral variations, so that a user, typically a geophysicist, 
can decide whether such frames should be discarded 
or retained for further processing. 

The present invention will be described further, 
by way of example only, with reference to an embodi- 55 
ment thereof as illustrated in the accompanying draw- 
ings, in which: 

Figure 1 shows a typical arrangement used to 



gather seismic data; 

Figure 2 is a diagram illustrating two different 
techniques for embodying the present invention; 
Figure 3 is a block diagram showing the elements 
of the apparatus according to the preferred em- 
bodiment of the present invention; 
Figure 4 is a flow diagram illustrating the process 
steps involved in dividing and re-ordering the 
edge list; 

Figure 5 is a diagram illustrating the re-ordering 
process carried out by the divider, and 
Figure 6 shows a frame of raw seismic data ot>- 
tained from a cable of sensors, and various other 
images representing various stages during the 
error detection process. 

Figure 1 illustrates a typical technique used to 
obtain seismic data representing the geological struc- 
ture of an area of rock below the earth's surface. 
Lines of cables 10 are set out on the earth's surface, 
each cable incorporating a number of geophone sen- 
sors 20. Then a source of sound 30 is placed within 
the vicinity of the cables and is activated to emit a 
shot, causing shock waves to pass down into the 
ground. These waves reflect off the various layers of 
rock 50 in the ground 40 and the reflected waves are 
detected by the array of geophone sensors 20. Next 
the source 30 is moved to a different position and 
again a shot is fired. The process is repeated a num- 
ber of times so as to create a series of shots at dif- 
ferent locations, and to yield a large number of frames 
of seismic data. 

This data is then returned to the laboratory for 
processing. For each cable the data received by the 
various geophone sensors on that cable for a partic- 
ular shot can be used to create an image 610 such as 
that shown in Figure 6. In image 610 the x-axis cor- 
responds to geophone position (ie trace) and the y- 
axis corresponds to time. Figure 6 shows a frame of 
raw seismic data obtained from a cable of sensors, 
and various other images representing various sta- 
ges during the error detection process; this figure is 
purely illustrative, and in the preferred embodiment of 
the present invention there is no need to represent the 
data graphically. Although figure 6 illustrates a rather 
idealised sample there will generally be an axis of 
symmetry similar to that shown. 

In computer vision literature, there are various 
known techniques that attempt symmetry detection, 
see for example the article entitled "Using a Mixed 
Wave/Diffusion Process to Elicit the Symmetry Set" 
by Scott, Turner and Zissermann, Proceedings of the 
Fourth Alvey Vision Conference, University of Man- 
chester, 31 August - 2 September 1988, page 221, 
and the article entitled "Recognize the Similarity Be- 
tween Shapes under Affine Transformation" by Hong 
and Tan, IEEE Second International Conference on 
Computer Vision, 5-8 December 1 988, page 489. Nei- 
ther of these techniques are particularly robust or ef- 
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f icient Applications in the computer vision field have 
generally been to simplified, artificial images. Seis- 
mic data presents a problem which is quite different 
in character since it generally has significantly great- 
er information content but does have the simplifies- s 
tion that the whole image exhibits symmetry about a 
fairly obvious axis that is usually vertical (since the 
rock interfaces tend to be approximately horizontal). 
Most known symmetry detectors employ complex 
techniques which concentrate on the problem of find- 10 
ing the axis of symmetry; clearly these techniques 
are not required for typical frames of seismic data. 

In the preferred embodiment of the present inven- 
tion the problem of finding the axis of symmetry can 
be solved either by using a straightforward peak de- is 
tector or by using information about the positioning of 
the shot with relation to the geophones. This will be 
discussed in more detail later with reference to Figure 
3. With the axis identified, the matching problem then 
becomes one of comparing the set of data to the left 20 
of the axis with the set of data to the right of the axis. 

Figure 2 is a block diagram illustrating at a high 
level two different techniques of implementing the in- 
vention to detect errors in a frame of seismic data. In 
both techniques the first step is to identify the axis of 25 
symmetry in the frame of seismic data represented by 
image 100. 

Then if technique (a) is chosen, the seismic data 
is divided into two sets of data; the first set contains 
data that can be represented by the left half of the i m- 30 
age 100, and the second set contains data which can 
be represented by the right half of image 100. The 
data in the first set is then reordered so that it would 
be represented by an image which is the mirror image 
of the image previously representing that set. The left 35 
set can then be represented by image 120 whilst the 
second set can be represented by image 130. 

These two sets of data are then passed separate- 
ly through an edge detector, typically employing 'Can- 
ny* edge detection, to yield two edge lists 140 and 40 
150. The edges in this context can be defined as dis- 
continuities in intensity which are associated with 
boundaries between areas in the original data or im- 
age representing that data. 

Finally the two edge lists are passed to a match- 45 
ing device 160 for comparison with each other. The 
matching device must then apply some kind of sym- 
metry matching method to the edge lists. Since one 
of the sets of data has been re-ordered (or 'mirrored'), 
the comparison that needs to be made can be seen 50 
to be similar to the matching of a pair of images in 
stereoscopic vision, where two cameras (or two eyes) 
see two slightly different views. There are some ro- 
bust algorithms available for stereo matching; the 
PMF algorithm produced by Pollard et al has shown 55 
itself to be particularly suitable. Details of this algo- 
rithm can be found in the article by S B Pollard, J E 
W Mayhew, and J P Frisby entitled "PMF: A Stereo 



Correspondence Algorithm using a Disparity Gradient 
Constraint", Perception, 14:449-470 (1985). 

If technique (b) is chosen the edge detection 
stage is carried out directly after the identification of 
the axis of symmetry, and before the frame of seismic 
data is divided. This yields an edge list 170 containing 
all the discontinuities located in the whole frame of 
seismic data. This edge list is then split using details 
of the axis of symmetry to produce two edge lists. The 
first edge list is then reordered using a process de- 
scribed later to yield a modified first edge list which, 
if represented graphically, would be a 'mirror im- 
age' of the representation of the unmodified first 
edge list These two edge lists 180 and 190 are then 
passed to the matching device 160 for comparison. 

Technique (b) is preferred over technique (a) 
since it reduces computing overhead by requiring 
only one set of data to be passed through the edge de- 
tection stage. 

The apparatus of the preferred embodiment of 
the present invention will now be described with ref- 
erence to the block diagram shown in Figure 3. 

The results from a particular set of seismic tests 
are stored in an external storage device 210 such as 
a magnetic tape device. If there is a general under- 
standing of the rock distribution in the survey area 
from which these results were obtained then some 
pre-processing correction can be made to the raw 
data. Since sound travels faster through matter as the 
matter's density increases, knowledge of the rock dis- 
tribution can be used to establish an understanding of 
the velocity distribution across the survey area, and 
this can be used to introduce a correction factor into 
the raw data stored in storage device 210. Once this 
correction has been made, images representing such 
seismic data will be expected to exhibit an inherent 
symmetry about an axis. This correction technique is 
one that is generally applied to such seismic data and 
so will not be described in any further detail herein. 

Once any pre-processing corrections have been 
made, then the apparatus 200 is activated to read one 
frame of seismic data from the external storage de- 
vice 210 into internal storage 220. This data will also 
include information about the width and depth of the 
data, and about the positioning of the shot with re- 
spect to the geophone sensors. 

The axis detector 240 then accesses the seismic 
data in order to select an axis as the axis of symmetry. 
In one embodiment the axis detector will employ a 
straightforward peak detection technique to deter- 
mine the position in the data which has the highest in- 
tensity; this position is then used as the x coordinate 
for a vertical axis of symmetry. 

However in the preferred embodiment the axis 
detector operates by locating the information in the 
data giving the positioning of the shot for that frame 
of data, and then uses that location as the x-coordin- 
ate to create a vertical axis of symmetry. 
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Next the frame of seismic data is passed to the 
edge detector 250, wherein a Canny edge detection 
routine is applied to the data to identify the data 
points representing discontinuities in the data (graph- 
ically, these are the lines of maximum gradient or 
most rapidly changing intensity). These data points 
are stored as an edge list 170, which is then passed 
to the divider 260; the divider also receives details of 
the axis of symmetry from the axis detector 240. 

The divider 260 divides the edge list into two sep- 
arate lists with data points to one side of the axis of 
symmetry being stored in the first edge list and data 
to the other side of the axis being stored in the second 
edge list Once this has been done a mirroring means 
within the divider 260 is then activated to reorder the 
first list The reordering is such that the reordered first 
edge list can be represented by an image which is the 
reflection (or •mirror') of the image which would rep- 
resent the unaltered first edge list The technique 
used to make this separation and reordering will be 
described in more detail later with reference to Figure 
4. 

The second edge list 190 and the reordered first 
edge list 180 are then passed to the matching device 
160, wherein the PMF stereo matching algorithm is 
applied to the edge lists. This yields statistics about 
the proportion of edge points in one edge list that are 
matched or unmatched to edge points in the other 
edge list In the preferred embodiment the matching 
algorithm also yields values for the average disparity 
and the average absolute disparity; for each pair of 
matched points, the disparity is the displacement be- 
tween the point in one edge list and its mirror point 
projected onto it from the other edge list. 

Then in the preferred embodiment the axis de- 
tector alters the x-coordinate of the axis of symmetry 
by one unit to the left (eg by the distance between two 
adjacent geophone sensors). The divider and match- 
ing device then repeat the above described procedure 
using the new axis of symmetry to yield a new set of 
statistics. The whole procedure is then repeated once 
more but with the axis of symmetry shifted by one unit 
to the right of the original axis of symmetry. 

Having done this the matching device will have 
created three average absolute disparity values, one 
for each axis of symmetry. If 'A was the original x-co- 
ordinate of the axis of symmetry and 'd* is the dispar- 
ity value then the matching device will have produced 
the following information: 



AXIS A 



AXIS A+l 
<*3 



If d, is numerically less than d 2 and d 3 , AXISA is 
the axis with the local minimum disparity value and 
the statistics generated by the matching device 160 



for AXISA are then passed on to decision logic 270. If 
however d 2 or d 3 have values less than di then an- 
other iteration of the process is carried out with a fur- 
ther axis of symmetry chosen one unit to the left or 

5 right of the axis with the lowest disparity value, eg if 
d 2 has a disparity value less than di and d 3 then an- 
other iteration is carried out with the axis of symmetry 
having an x-coordinate of f A-2\ If this yields a dispar- 
ity value greater than d 2 then d 2 is the local minimum 

10 value and the statistics generated by the output from 
the matching device 160 when the axis of symmetry 
had an x-coordinate of A-1 are passed to the decision 
logic 270. 

Once the axis of symmetry giving rise to the local 

is minimum disparity value has been identified, the de- 
cision logic 270 will determine whether the statistics 
(eg disparity value, percentage of edges matched, 
etc) output by the matching device are within thresh- 
olds preset by the user. If they are then the frame of 

20 data stored in internal storage 220 can be passed di- 
rectly to the data processing stages 280 since the in- 
dication is that the frame of data is reasonably free 
from errors and hence suitable for processing. If how- 
ever the statistics are greater than the threshold val- 

25 ues then the frame of data is referred for user eval- 
uation 290 since the indication is that it includes an 
unacceptable proportion of errors. The user will typi- 
cally be a skilled geophysictst and will view an image 
representing the frame of data to see 

30 how 'good* or 'bad' it looks. If the user agrees that the 
data is 'bad' then he can discard that frame of data. 
If on the other hand he decides that the data looks ac- 
ceptable he can pass the frame of data on to the data 
processing stages 280. 

35 Ideal threshold values for the statistics are typi- 

cally set by the user during an initial tuning phase. In 
this initial tuning phase the user will view a number of 
frames of seismic data using the standard subjective 
prior art approach and decide whether each frame 

40 is 'good 1 or 'bad'. These same frames of data will 
then be passed one by one through the apparatus of 
the preferred embodiment to see if the decision logic 
270 refers any of them for user evaluation. Ideally 
only the ones considered by the user to be bad should 
45 be referred to him, but typically the threshold values 
will have to be manipulated until some compromise is 
struck (eg all of the bad frames are passed for user 
evaluation plus a small proportion of the good 
frames). 

so The decision logic can be adapted to take ac- 

count of a number of features when making its deci- 
sion as to the quality of the frames of data. It has been 
found that the disparity values created by the PMF al- 
gorithm are good statistics to use as part of this deci- 

55 sion process. If the average absolute disparity fails to 
reach a minimum in the region of the expected axis 
of symmetry then it is likely that an error exists; this 
fact can be used by the decision logic 270. In addition, 
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for good data, the minimum disparity value and the 
depth of the average absolute disparity curve should 
fall below chosen threshold values. Further the dis- 
parity curve should not exhibit ancillary minima. Giv- 
en this information the decision logic has shown to s 
perform well. 

By using the above described apparatus many 
images can be inspected without manual interven- 
tion. If the thresholds have been chosen correctly 
very few images will be referred for manual inspec- 10 
tion. The user can decide the optimal balance for the 
threshold values in order to minimise missed errors 
while maximising throughput 

The above technique in which vertical axes are 
chosen as candidate axes of symmetry seems to is 
work well when the seismic data represents reflec- 
tions off rock interfaces lying roughly horizontal. If 
however the rock interfaces were at a significant an- 
gle to the horizontal then frames representing such 
interfaces may well not give good matching and would 20 
be referred to the user. The user would then notice 
that the data referred to him was in fact 'good' data 
but that the axis of symmetry was slightly skewed 
from the vertical, and would hence pass the data on 
for subsequent data processing rather than discard- 25 
ing it 

Alternatively some form of skew correction circui- 
try could be developed to apply a correction at an ap- 
propriate stage in the process. 

The process of dividing and reordering the edge 30 
list will now be described in more detail with reference 
to the flow diagram shown in Figure 4. 

At step 41 0 an axis of symmetry is created using 
information contained in the frame of seismic data 
about the positioning of the sound source during that 35 
shot; this process is carried out by the axis detector 
240 described earlier with reference to Figure 3. 

At step 420 an edge list is created from the frame 
of seismic data by performing Canny edge detection 
on that frame. The output of the Canny edge detection 40 
is an edge list listing all of the data points at which dis- 
continuities occur. 

Then at step 430 the first data point is accessed 
from the list and its x-coordinate value is compared 
with the x-coordinate of the axis of symmetry as de- 45 
termined at step 410. In the preferred embodiment 
the comparison involves subtracting the x-coordinate 
of the axis of symmetry from the x-coordinate of the 
data point If the result of the subtraction is positive 
(step 440) the data point is written directly to the sec- 50 
ond edge list (step 450). The positive value indicates 
that the edge point lies to the right of the axis of sym- 
metry, and the value resulting from the subtraction 
gives the data point's position as adjusted to reflect 
its offset from the axis of symmetry. 55 

If however the result is negative the process pro- 
ceeds to step 460 where it is determined whether the 
current data point is on a different row to the previous 



data point processed. If it is not then the current data 
point is written to a temporary buffer (step 470). If 
however the current data point is on a new row then 
the contents already in the temporary buffer are read 
out in reverse order and stored in the first edge list 
(step 480). The buffer is then cleared (step 490) and 
the current data point then written to the buffer at step 
470. This ensures that only the contents of one row 
are stored in the buffer at any one time. 

The reason for writing such data points (for which 
the result of the subtraction is negative) to a buffer 
rather than directly into the first edge list is as follows. 
The negative value indicates that the data point lies 
to the left of the axis of symmetry as viewed in Figure 
6. It is preferable to perform a reordering (or •mirror- 
ing") operation on all data points to the left of the axis 
of symmetry in order to reduce the complexity of later 
comparison of the left-side edge list with the right- 
side edge-list. However, for each left-side edge point 
this involves modifying the edge orientation as well 
as adjusting the position with respect to the axis. Sim- 
ple "mirroring 11 by itself is insufficient as the compar- 
ison process relies to some extent on the ordering of 
the edge-list Accordingly, for the left-side edges, the 
edge points must be re-ordered within each row of the 
data. To do this it is necessary to collect all of the left- 
side edge points for any given row together in a buffer 
and then to perform the reordering. 

At step 500 it is determined whether all of the data 
points in the edge list have been compared with the 
axis of symmetry. If this is not the case then the proc- 
ess returns to step 430 and a further subtraction step 
is performed on one of the remaining edge points. If 
all of the edge points have been compared, and so 
have been written either to the second edge list first 
edge list, or the temporary buffer, then the process 
proceeds to step 510. At this step the data points re- 
maining in the buffer (corresponding to the final row) 
are read out of the buffer in reverse order and stored 
in the first edge list The process has then achieved 
the necessary reordering discussed above. Finally 
the buffer is cleared (step 520) and both edge lists are 
passed to the matching device 160 (step 530). 

Figure 5 is a diagram illustrating this re-ordering 
process. The frame of seismic data is represented by 
image 550, which shows the nth time slice (or row) in 
detail. A portion of the edge list for the nth row is rep- 
resented by chart 560. Here we can see that four edge 
points C1, C2, C3, C4 have been identified by the 
edge detector for row n. 

The subtraction step as illustrated in figure 5 
shows that subtraction of the x-coordinate of the axis 
of symmetry from the x-coordinate of the edge points 
of row n yields negative values for C1 and C2 and pos- 
itive values for C3 and C4. The positive values are 
written directly into the second (or right-side) edge 
list a portion of which is shown by chart 580, whilst 
the modulus of the negative values are written into 
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the buffer 570. 

The re-ordering step involves, upon reaching the 
first left-side edge point in row n+1, reading the val- 
ues for row n out of the buffer 570 in reverse order and 
storing them in the first (or left-side) edge list, a por- s 
tion of which is shown by chart 590. In this way the 
first edge list comes to be filled with data points which 
can be represented by an image which is the reflec- 
tion (or 'mirror') of the image which would represent 
a first edge list created without any re-ordering. 10 

Figure 6 shows a frame of raw seismic data ob- 
tained from a cable of sensors, and various other im- 
ages representing various stages during the error de- 
tection process. As mentioned earlier these images 
are purely illustrative and do not need to be repro- 15 
duced during typical operation of the apparatus of the 
preferred embodiment 

Image 610 shows a frame of raw data obtained 
from one cable of geophone sensors for one particu- 
lar shot; this image will preferably have been correct- 20 
ed for noise using a known noise correction techni- 
que. Image 620 is illustrative of the image that can be 
produced to represent the edge list created by edge 
detector 250. The axis of symmetry as generated by 
the axis detector 240 is shown superimposed on im- 25 
age 620. 

As discussed the divider takes the edge list rep- 
resented by image 620 and creates right and left-side 
edge lists, the data in the left-side list having been re- 
ordered; these two edge lists are represented by im- 30 
ages 630 and 640 respectively. Image 650 illustrates 
the proportion of data points matched by the match- 
ing device 160, whilst image 660 illustrates the pro- 
portion of unmatched data points. 

The value of the above described preferred em- 35 
bodiment (corresponding to technique (b) in Figure 2) 
is that it reduces computing overhead by performing 
a single edge-detection step, followed by an efficient 
splitting process. Its effectiveness has been demon- 
strated by processing a seismic shot which was con- 40 
trived to be symmetrical and achieving 100% match- 
ing between the two sides. 



Claims 43 

1. Apparatus for detecting errors in frames of seis- 
mic data that can be presented as images, such 
images being expected to exhibit an inherent 
symmetry about an axis, the apparatus compris- so 
ing: 

a storage device for storing a frame of seismic 

data which can be presented as an image; 

a selection device for selecting an axis within the 

image as the axis of symmetry; 55 

a divider for dividing the data about the axis of 

symmetry to create first and second sets of data; 

and 



a comparison means for comparing the sets of 
data according to preselected matching criteria to 
determine, based on the degree of matching, 
whether the extent of mismatch in the frame of 
data is within a predetermined tolerance. 

2. Apparatus as claimed in claim 1, wherein the 
comparison means comprises: 

a matching device for comparing the sets of data 
to determine the degree of matching between 
such sets; and 

decision logic, responsive to output from the 
matching device, to determine whether the extent 
of mismatch is within the predetermined toler- 
ance. 

3. Apparatus as claimed in claim 2, wherein the 
comparison means further comprises: 

an edge detector for determining the data points 
representing discontinuities in the image, and to 
create an edge list of such data points, such edge 
lists being used as the basis of comparison by the 
matching device. 

4. Apparatus as claimed in claim 3, wherein the div- 
ider operates on the whole frame of seismic data 
prior to passing the first and second sets of data 
separately through the edge detector, the edge 
lists created by the edge detector being received 
by the matching device for comparison. 

5. Apparatus as claimed in claim 3, wherein the 
edge detector operates on the whole frame of 
seismic data and then the edge list is operated on 
by the divider so as to create two edge lists as the 
first and second sets of data for comparison by 
the matching device. 

6. Apparatus as claimed in any preceding claim, 
wherein the divider includes a mirroring means 
for re-ordering the data in the first set of data so 
that it can be represented by an image which 
would be the mirror image of the image previous- 
ly representing that set, the comparison means 
then comparing the second set of data with the 
re-ordered first set of data. 

7. Apparatus as claimed in claim 5 and claim 6, 
wherein the mirroring means re-orders the first 
set of data created by the divider by employing a 
buffer in which to store data from that set of data 
constituting a row of an image, the contents in 
the buffer then being read out in reverse order 
and stored as the re-ordered first set of data. 

8. Apparatus as claimed in any preceding claim, 
wherein the selection device is an amplitude 
peak detector wh ich select s the axis of symmetry 
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as the position in the data having the highest de- 
tected intensity. 

9. Apparatus as claimed in any of claims 1 to 7, 
wherein the selection device uses information s 
contained in the seismic data, giving the position- 
ing of the shot for that frame of data, to select an 
axis as the axis of symmetry. 

10. Apparatus as claimed in any of claims 3 to 9, 10 
wherein the edge detector employs a canny edge 
detection algorithm in order to identify the data 
points representing edges. 

11. Apparatus as claimed in any preceding claim, is 
wherein the matching device uses a stereo 
matching algorithm. 

12. Apparatus as claimed in claim 11, wherein the 
stereo matching algorithm is the PMF algorithm. 20 

13. Amethod of detecting errors in frames of seismic 
data that can be presented as images, such im- 
ages being expected to exhibit an inherent sym- 
metry about an axis, the method comprising the 25 
steps of: 

(a) storing a frame of seismic data that can be 
presented as an image in a storage device; 

(b) selecting an axis within the image as the 

axis of symmetry; 30 

(c) employing a divider to divide the image 
about the axis of symmetry to create first and 
second sets of data; and 

(d) comparing, in a comparison means, the 
sets of data according to preselected match- 35 
ing criteria to determine, based on the degree 

of matching, whether the extent of mismatch 
in the frame of data is within a predetermined 
tolerance. 

40 

14. A method as claimed in claim 13, wherein the 
comparison step (d) comprises the steps of. 

(i) comparing, in a matching device, edge lists 
to determine the degree of matching between 

the edges represented by such lists; and 45 

(ii) determining, in response to output from the 
matching device, whether the extent of mis- 
match is within the predetermined tolerance. 

15. A method as claimed in claim 14, wherein the so 
comparison step (d) further comprises, prior to 

the comparing step (i), the step of: 
(Hi) determining, in an edge detector, the data 
points representing discontinuities in the image, 
and creating an edge list of such data points, ss 
such edge lists being used as the basis of com- 
parison by the matching device. 



16. A method as claimed in any of claims 12 to 15, 
wherein the dividing step (c) is performed prior to 
any of the steps in comparison step (d). 

17. A method as claimed in any of claim 15, wherein 
the edge detection step (d)(iii) is performed on 
the whole frame of data and then the edge list is 
used as the input for step (c) so as to create two 
edge lists as the first and second sets of data for 
comparison at step (d)(i). 

18. A method as claimed in any of claims 13 to 17, 
wherein the dividing step (c) includes a mirroring 
step which re-orders the data in the first set of 
data so that it can be represented by an image 
which would be the mirror image of the image 
previously representing that set, the comparison 
step (d) then comparing the second set of data 
with the re-ordered first set of data. 
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(g) Error detection In seismic data, 

(57) The present invention relates to the detection 
of errors in seismic data and in particular to 
apparatus for detecting errors in frames of seis- 
mic data that can be presented as images. Due 
to the techniques used for collecting seismic 
data, an image produced to represent that data 
will typically be expected to exhibit an inherent 
symmetry about an axis. The apparatus of the 
present invention comprises a storage device 
for storing a frame of seismic data which can be 
presented as an image, a selection device for 
selecting an axis within the image as the axis of 
symmetry, and a divider for dividing the data 
about the axis of symmetry to create first and 
second sets of data. Further it includes a com- 
parison means for comparing the sets of data 
according to preselected matching criteria to 
determine, based on the degree of matching, 
whether the extent of mismatch in the frame of 
data is within a predetermined tolerance. 
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